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ABSTRACT 

IRAS 16279—4757 belongs to a group of post-AGB stars showing both PAH 
bands and crystalhne sihcates. We present mid-infrared images, that resolve 
the object for the first time. The morphology is similar to that of the 'Red 
Rectangle' (HD 44179), the prototype object with PAHs and crystalline sili- 
cates. A two-component model and images suggest a dense oxygen-rich torus, 
an inner, low-density carbon-rich region and a carbon-rich bipolar outflow. The 
PAH bands are enhanced at the outflow, while the continuum emission is con- 
centrated towards the center. Our findings support the suggestion that mixed 
chemistry and morphology are closely related. We discuss the ISO/SWS spectra 
of IRAS 16279— 4757. Several bands in the ISO/SWS spectrum show a match 
with anorthite: this would be the first detection of this mineral outside the solar 
system. Compared to HD 44179, the shapes of PAH bands are closer to those 
of planetary nebulae, possibly related to a population of small PAHs present 
HD 44179, but absent around IRAS 16279— 4757. Detailed examination of the 
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spectra shows the individual character of these two objects. The comparison 
suggests that the torus found in IRAS 16279—4757 may have formed more re- 
cently than that in HD 44179. 

Subject headings: stars: AGB and post-AGB, circumstellar matter — infrared: 
stars 



1. Introduction 

Some post-AGB stars show both PAH and crystaUine sihcate bands in their infrared 
spectra. The formation history of this mixed chemistry (oxygen-rich silicates versus carbon- 
rich PAHs) is not well understood. A possibihty, but implausible, is that these stars all 
evolved from oxygen-rich to carbon-rich within the last few hundred years (Zijlstra ct al. 
1991). Waters ct al. (1998) and Molstcr ct al. (1999) propose that the silicate dust is stored 
in a long-lived circumbinary disk. In this scenario, the PAHs form during a later mass- 
loss phase, after the star became carbon-rich, while the gas stored in the disk retains the 
chemistry of the earlier, oxygen-rich phase. Part of the amorphous silicate dust crystallizes 
in the disk. The scenario explains why there are relatively few post-AGB stars with this 
mixed chemistry, and why the silicates have a lower temperature than the carbon rich dust. 
It requires all such stars to be binaries, as is the case for the prototype of the class, the Red 
Rectangle (HD 44179; Waelkens et al. (1996)). 

IRAS 16279-4757 (hereafter IRAS 16279) is classified as a post-AGB star, based on a 
double peak in its spectral energy distribution indicating a detached envelope (van der Veen 
et al. 1989). Optical spectra suggest a spectral type of G5 (Hu et al. 1993). The ISO/SWS 

spectra show crystalline silicates beyond 20 (Molster et al. 1999) and PAH bands are seen 
in the near-infrared (van der Veen et al. 1989) and ISO spectra. IRAS 16279 is therefore a 
member of the group of mixed-chemistry post-AGB stars. 

We present TIMMI-2 mid-infrared imaging and spectroscopic data, resolving IRAS 16279 
for the first time. In this paper we discuss the spatial distribution of the different dust compo- 
nents based on these images, and compare these with predictions from the circumbinary disk 



^ Based on observations with the European Southern Observatory, 3.6-meter telescope with TIMMI-2 at 
La Silla. The proposal number is 71.D-0049. 

^Based on observations with ISO, an ESA project with instruments funded by ESA Member States (espe- 
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of ISAS and NASA. 
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scenario. We compare the spectra of this object with those of the prototype mixed-chemistry 
object, the Red Rectangle. 

2. The observations 

We obtained 10 and 18 //m images with TIMMI-2 (Kaufl et al. 2000) at the ESO 3.6m 
telescope on La Silla, on the 27th of May 2003. The weather was clear with occasional cirrus. 
We used 6 narrow N and Q-band filters (Table 1). The diffraction hmit is about O.Tarcsec 
(first null) in the N-band and 1.4arcscc at Q-band. The optical seeing of 0.8-1.1 arcsec in 
FWHM allowed us to reach the diffraction limit in the images. The pixel scale is 0.2 arcsec. 
The background was subtracted with nodding and chopping with an offset of 15 arcsec. The 
pipeline-reduced data were used for imaging. The results from the aperture photometry are 
given in Table 1. The exposure time is about 16 minutes per image. 

We also obtained 10 fim spectra with TlMMl-2, taken through a 1.2x50 arcsec^ slit, 
oriented North-South. The slit was centered 0.5 arcsec west from the brightest N-band 
position. The resolution is about 160 in A/ A A. The background was subtracted with 15 
arcsec chop-and-nod. We used eclipse and IDL for the data reduction. The source extends 
approximately 6 x 8 arcsec^, and the spectra are spatially resolved. The pixel scale of 0.45 
arcsec per pixel was binned to 0.9 arcsec along the slit to increase the signal-to-noise ratio. 
The total exposure time is 58 minutes. 

The flux calibration, and the response/atmospheric transmission correction in the spec- 
tra was derived using HD 169916 (KlIII). The spectral template for HD 169916 is calculated 
by Cohen (1998). 

Table 1: Filters and flux densities. The last column lists the aperture used for photometry. 



Filter 


Ar.o'Xcui [/'111] 


Flux [Jy] 


Error [J}-] 


Apt . [arcsec] 


N7.9 


7.42- 8.11 


23.4 


4.8 


4 


N8.9 


8.29- 9.07 


19.5 


4.4 


4 


N9.8 


9.10-10.02 


20.0 


4.5 


4 


N10.4 


9.80-10.82 


29.0 


5.4 


4 


Nil. 9 


10.99-12.19 


59.8 


7.7 


6 


Q 


17.35-18.15 


172. 


12. 


6 
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We also use archival ISO/SWS spectra, obtained on 21 August 1997 with a resolution 
of A/AA =250-600. The aperture is 14x27arcsec^ or larger, sufficient to cover the entire 
source. However, there are several flux discrepancies between the detector bands, probably 
due to a slight mis-pointing. We scaled the separate ISO/SWS bands (12.45-27.29 /um, 
27.28-28.90 fim, and 28.90-45.38 fim) so as to smoothly connect across the bands, and to 
simultaneously agree with the IRAS broad-band flux levels. 

3. Results 

3.1. Mid-infrcired images 

The TIMMI-2 images of IRAS 16279 are shown in Fig. 1. The object is clearly resolved 
in all images, and more elongated within the PAH bands. This is clearly seen in Fig. 2 
in which we compare the intensity profiles of the various filters along the north-south axis 
with the profiles of a calibration star. To further quantify the spatial variations in the PAH 
and continuum emission, we show the PAH-to-continuum-ratio (Figs. 3 and 4). The PAH 
emission is enhanced towards the SW and NE, giving the appearance of bipolarity. At the 
center, the PAH-to-continuum ratio is low. The elongated PAH emission distribution, and 
the low PAH / continuum ratio at the center are similar to what is found in the Red Rectangle 
(Waters et al. 1996). 

To confirm that the PAH images are more elongated than the continuum images, we 
fitted the images with two 2-dimensional Gaussian functions. The first Gaussian component 
represents the bright core, and the second Gaussian measures the intensity and elongation 
direction of the faint, extended emission. The peak positions of the two Gaussian profiles 
are fixed at the same place. The best- fit parameters are summarized in Table 2. The PAH 
bands (N7.9, N8.9, Nil. 9) show a relatively stronger second component (larger 02/01 ratio), 
confirming the extended appearance in these bands. 

The Q-band and Nil. 9 images (Fig. 1) show a more rectangular shape, slightly elongated 
both along the SW-NE direction (as found in PAH images) and the SE-NW direction. The 
rotation angles change from 160 degrees in the 10-//m bands, to 75 degrees at 11.9 /xm, to 10 
degrees at 18 //m. The Q-band filter contains mainly dust continuum (Sec. 3.3). Although 
the Nil. 9 band contains a PAH feature, this contributes only about 1/2-1/3 of the intensity, 
and the rest is continuum. This rectangular shape might also be present in other continuum 
images, however, the emission is fainter in the shorter-wavelength continuum bands. 
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Fig. 1. — The images of IRAS 16279—4757 in units of logF,^ in Jy per pixel. The object is 
elongated, especially in the PAH bands. Contour lines show 0.0025, 0.005, 0.01, 0.05, and 
O.lJy per pixel for the N-band, and 0.01, 0.025, 0.05, 0.1, 0.25, and 0.5 Jy per pixel for the 
Q-band. The dashed lines in (e) show the approximate slit position for the spectra. 



3.2. TIMMI-2 Spectra 

The N-band spectra (Fig. 5) show PAH bands at 7.7, 8.6, and 11.2 /im. The weaker 
12.7 yum PAH band (e.g. Cohen, Tielens, AUamandola 1985; Hony et al. 2001) is also detected 
in IRAS 16279. The spatially resolved N-band spectra show that the central spectrum (ac- 
tually 0.5" west of the brightest position) has the strongest continuum, while the strongest 
11.2 /im PAH band is found 0.9" south (see also Fig. 6). The spectrum 0.9" north shows only 
a weak PAH feature. The line-to-continuum ratio is lower at the center than at 0.9" south, 
in agreement with the ratio images displayed in Fig. 3. 

The lower PAH-to-continuum ratio at the center (Figs. 3 and 4) is due to the fact that 
the continuum intensity declines sharper than the PAH intensity, along the SW-NE direction 
(Fig. 6). The PAH intensity itself is still reasonably high near the center as seen in spectra. 
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3.3. Modeling the ISO/SWS spectra 

In order to obtain a globally consistent physical picture of the source we model the 
ISO/SWS spectrum. We show the full 2.3-45 ;um spectrum in Fig. 7. Most of the detected 
energy falls in the infrared range and thus originates from the circumstellar dust. The peak 
wavelength of the energy distribution (~30 /xm) implies a typical dust temperature of about 
100 K. Below 15 //m the spectrum is dominate by C-rich material, while beyond 20 //m the 
solid-state features can be attributed to crystaUine silicates. 

We use a spherically symmetric model consisting of two distinct n oc regions, one 
containing carbon dust and PAHs, and the other containing silicate dust. The model is 
described in Siebenmorgen, Kriigel (1992), and Siebenmorgen, Zijlstra, Kriigel (1994): it is 
a 1-d model which includes scattering and transient heating. The radiative transfer is solved 
independently for the carbon-rich and the oxygen-rich region. We do not include the possible 
effect of the fairly low extinction within the carbon-rich region on the heating of the (outer) 
oxygen-rich region. 

The distance and extinction are derived assuming a luminosity of 1 x 10^ L©. The 
integrated flux (optical to infrared photometry from Hu et al. (1993), Van de Steene et al. 
(2000), IRAS, and ISO) gives a distance estimate of 2.0 kpc. IRAS 16279 is in the Galactic 
plane {b = +0?09) and will suffer interstellar extinction. The optical photometry is reported 
as i? = 18.42 and I = 14.57 mag (Hu et al. 1993). For a spectral type of G5 (Hu et al.) 
we derive Ay = 14.8 mag. Six field stars within 30arcmin show an interstellar extinction of 
1.0 mag kpc~^ out to 1.3 kpc. This suggests an interstellar extinction of Ay ~2.0mag. Most 
of the large extinction of IRAS 16279 is likely circumstellar origin. Therefore, we ignore the 
interstellar extinction in the model fit. 

Fig. 7 shows the fitted result. Up to 10 /xm, the emission from carbon material is the 
dominant source of continuum emission. The oxygen shell begins to contribute at 11 /im and 
is dominant beyond 15 //m. The derived parameters (Table 3) show that most of the dust is 
actually cold oxygen-rich dust: the carbon region is much less dense and contributes little 
to the mass. 

Table 3 shows that the carbon-rich gas is mostly located inside the oxygen-rich region. 
Both regions have a covering factor (sky coverage as seen form the central star) of about 30- 
35 %. The precise value depends on the choice of parameters, but all fits we obtained required 
that the output infrared flux be scaled down by at least half. The derived covering factors 
and the bipolar morphology in the images suggest that the different dust components may 
be located in different directions from the star. The model is consistent with an oxygen-rich 
torus, where we assume that the radial density distribution is the same everywhere within 
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the torus. 

The major PAH bands can be reasonably fitted with the model, except for the one at 
8.6 iim (this problem is also found by Siebenmorgen et al. (1994)). The oxygen shell gives less 
than half the observed Ay and a weaker-than-observed amorphous silicate absorption feature 
(seen against the carbon-continuum): the remainder may be due to interstellar absorption, 
or may be due to geometrical effects. The 10 /xm absorption is seen mainly against the inner 
arcsecond, carbon-rich dust. For a shell, the tangential line of sight can have a much 
higher extinction, and depending on the viewing angle, this can affect the carbon-rich region. 
The 1-d nature of the model means that self-absorption within the dust regions is not well 
described. The line of sight to the star also shows higher extinction, but it is clear that a 
global 1-d model may well underestimate the column density in individual lines of sight. 

The carbon rich shell has an outer radius of 5 x 10^^ cm (1.8 arcsec) in the model, which is 
close to (T2maj measured in the PAH bands. The inner radius of carbon-rich shell is 1 x 10^^ cm 
(0.3 arcsec). The inner radius is not clearly seen in the image because the spatial resolution 
is not high enough and because the point spread functions of the central star and the shell 
are merged. 

The model fitting yields several constraints: 

(1) To excite the PAH, the model requires a stellar temperature of 6 000 K or more; other- 
wise the bands disappear due to lack of exiting photons. This includes absorption by 
neutral PAHs. The PAH emission features are thought to arise from mixture of neutral 
and charged PAHs. PAH ions can absorb in the optical (c.f. Li, Drainc 2002) and this 
would allow a lower stellar temperature. Li, Draine (2002) show that weak PAH bands 
can be seen for stellar temperatures as low as 3000 K, but reach full strength only 
above 6000 K. Additionally, the presence of 11.2 /xm feature, which attributes to neu- 
tral PAHs (Hony et al. 2001) suggests that UV radiation is necessary and extremely low 
temperature is not possible. The observed PAH bands are compatible with a G-type 
central star. 

(2) The PAH bands are only seen in relatively unshielded (low extinction) gas. Within 
the model, only a low density shell contributes to the PAH bands. Higher density 
carbon regions, if present, would not contribute to the observed PAH bands: the 
excitation of the PAHs would be suppressed by internal extinction. We therefore have 
few constraints on higher density carbon-rich gas. 

(3) The oxygen-rich gas is located relatively far out (inner radius of 1 arcsec) and has much 
higher mass than the carbon-rich gas. The low temperatures of this gas are consistent 
with the observed crystaUine sihcate bands (70-110 K; Sect. 4.2). 
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4. Dust features 

The N-band images clearly show that the PAH emission originates from a bipolar ge- 
ometry. This further strengthens the link between this source and HD 44179, and the link 
between the mixed chemistry nature and a bipolar (+disk) geometry. It is interesting to fur- 
ther explore the comparison between IRAS 16279 and the Red Rectangle as the prototype 
of the dual chemistry sources. In the following we will compare the PAH emission features 
between the two sources, and we present a detailed comparison of the crystalline silicate 
bands. 



4.1. PAH bands 

4.1.1. Integrated spectrum 

In Fig. 8 we show the ISO PAH spectra of IRAS 16279, HD 44179 as the prototype of the 
mixed chemistry post-AGB stars and CD-42 11721 as a prototypical ISM/star-forming- region 
PAH spectrum. The profile features are very similar in the 3 fim and the 10 to 13 fim region 
(C-H modes) while large differences are found in the PAH spectrum from 6 to 9 fim (C-C 
modes). This strong variation in the 6 to 9 /xm PAH bands when compared to the other PAH 
bands is a common feature of the interstellar PAH spectra (van Diedenhoven et al. 2004). 
The profiles of the 6-9 //m PAH bands of IRAS 16279 are intermediate between HD 44179 
and CD-42 11721 (Peeters et al. 2002). The 6.2 pim feature of IRAS 16279 is broader than 
in CD-42 11721. This broadening is due to an extra contribution near 6.05 //m which is also 
present in HD 44179 and a weak emission band which peaks near 6.28 fim. The latter band is 
dominant in the spectrum of HD 44179 while only marginally detected in IRAS 16729, where 
the usual 6.22 /xm band is the strongest. 

The PAH bands from 7 to 9 fira (regularly dubbed the "7.7 complex") also attest to 
the intermediate nature of IRAS 16279 PAH emitters. In the ISM and in star-forming re- 
gions, usually the 7.6 fim band is dominant, while in planetary nebulae and post-AGB stars 
the 7.8 /xm band is strongest (Bregman 1989; Cohen et al. 1989; Peeters et al. 2002). In 
IRAS 16279 both bands are roughly equally strong. 

Hony et al. (2001) found that the band strength ratio of the 11.2 to 12.7/xm bands 
(/ii,2//i2.7 where In. 2 and /12.7 represent the integrated flux of each band) correlates with 
the source type. In the planetary nebulae this ratio is very high as is the case for HD 44179 
(-^ii.2/A2.7=5.0). The band strength ratio found in IRAS 16279 (/ii.2/-?^i2. 7=2.0) is not typical 
for evolved stars. 
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Table 2: The fitted results of the images with two 2-dimensional Gaussian profiles, a is the 
peak intensity in Jy per pixel. (Tmaj and o"min are along the major and minor axes in arcsec, 
respectively. 6 is the rotation angle of the major axis with respect to North (clockwise). The 
subscripts 1 and 2 show each Gaussian component (ai > 02). The first component measures 
the bright core region, and the second component shows the more extended emission. The 
first component in N7.9 is strongly affected by the extended component, and 9i is totally 
different from those in other bands. 







chimin 




9i 


C"2maj 


C2mm 


a2 


02 


02 / ai 


N7.9 


1.28 


0.84 


0.122 


152.7 


1.79 


0.84 


0.034 


136.9 


0.279 


N8.9 


0.66 


0.40 


0.117 


79.5 


1.64 


1.06 


0.056 


158.0 


0.479 


N9.8 


0.73 


0.57 


0.163 


82.0 


1.84 


1.35 


0.025 


166.5 


0.153 


N10.4 


0.74 


0.58 


0.198 


81.9 


1.67 


1.29 


0.048 


167.3 


0.242 


N11.9 


0.69 


0.52 


0.322 


86.0 


1.64 


1.17 


0.135 


74.5 


0.419 


Q 


1.10 


0.79 


0.266 


146.1 


2.39 


1.73 


0.213 


8.4 


0.801 



Table 3: Model parameters of the two-shell model. 







C shell 


sheU 




[cm] 


1 X 10^6 


4 X 10^6 


^out 


[cm] 


5 X 10^^ 


8 X 10^^ 


-^shell 


[Mo] 


3 X 10-3 


2.1 


Av 


[mag] 


0.6 


6.4 


7d 


[K] 


330-140 


125-57 


Goverage 


[%] 


30 


35 
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When comparing the band strength of the 3.3 /xm band with the longer wavelength PAH 
bands we find that the 3.3 /xm of IRAS 16279 is much weaker than in HD 44179. The ratios 
of the integrated band strengths, /11.2//3.3, are 1.9, 5.2 and 3.9 in HD 44179, IRAS 16279, 
CD-42 11721 respectively. The ratios of IRAS 16279 are closer to CD-42 11721. This may 
point to a population of smaller PAHs around HD 44179 which does not exist around the 
other two sources. Such small PAHs might be easily destroyed when exposed to the hard 
UV radiation present in the ISM and planetary nebulae. 

4.1.2. Spatially resolved spectroscopy 

In Fig. 6 we show the profile of the 11.2 fim band as extracted from the N-band spectra 
for different positions. On the right we compare the normalized profiles. The profile varies 
with the position along the slit, being broader at the central position and 0.9" north. A 
longer tail at 11.5 //m is usually interpreted as being due to hotter PAHs (Tielens et al. 1999; 
Pech, Jobhn, Boissel 2002). Alternatively, a different profile might arise from a variation 
of the PAH population as a function of distance to the source. For example, PAHs with 
two neighboring C-H bonds may give rise to emission in the red wing of the 11.2 /xm solo 
C-H mode (Hony et al. 2001). We prefer the hot PAH interpretation, because of the fact 
that the sharp blue-side rise stays the same in the all spectra, which is difficult to envisage 
in the scenario of a strongly varying PAH population. Furthermore, the most pronounced 
tail is found closest to the central star, where multiple photon heating may play a role and 
therefore hotter PAHs might be expected. 

Miyata et al. (2001) show spatially resolved N band spectra of HD 44179 in which they 
also find variations of the emission features as a function of distance to the star. Like 
in IRAS 16279 the 11.2 /xm band diminishes close to the central star. However, unlike in 
IRAS 16279, the HD 44179 spectrum near the central star is dominated by bands at ~11.0 
and ~11.9//m. These bands could be related to the crystaUine enstatite that is present in 
the disk. Due to this possible confusion between the different contributors to the emission in 
HD 44179 it is not feasible to determine the band profile of the 11.2 //m nearest the central 
star. 

Concluding, we can state that there are significant differences in the PAH spectrum of 
HD 44179 and IRAS 16279, both in band strength ratios and profile features, where the latter 
source exhibits an PAH spectrum intermediate between HD 44179 and the prototypical ISM 
PAH spectrum. In this sense the PAH spectrum of IRAS 16279 compares better, in bands 
strengths and feature profiles, with the young C-rich planetary nebulae IRAS 21282+5050 
(Hony et al. 2001; Peeters et al. 2002). 
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4.2. Crystalline silicate bands 

The ISO/SWS spectra beyond 20 fim show the sohd state features which are attributed 
to crystalhne sihcates (Fig. 7). Continuum subtracted spectra are shown in Fig. 9. The 
continuum is obtained with a quadratic fit (polynomial fit to the three neighboring points) 
to the smoothed spectra. 

The feature at 33.7;um has been identified with crystalline olivine (Mg2a;Fe2-2a:Si04, with 
X between and 1). Its position fits best with that of the pure Mg-olivine, i.e. forsterite 
(Mg2Si04) (Waters et al. 1996; Molster et al. 2002a). The band ratio depends on the temper- 
ature (Molster 2000, chap. 1). The absence of a feature around 23.7 //m shows that forsterite 
is below 100 K. The plateau up to 37 /im, which has no clear mineralogical identification, is 
prominent in this star. 

The feature at 40.5 iim is normally attributed to pyroxene-like sihcates (MeSiOs, where 
Me stands for metal). The presence of features also around 43 /im is normally an indication 
for the presence of enstatite (MgSiOa). However, the lack of enstatite features below 30 /xm 
makes this identification questionable. The ISO/SWS band 3E (27.5-29 //m) is known for 
its mediocre performance. This could be the reason for the lack of enstatite features around 
28 fim. Alternatively, enstatite can also simply be very cold just like the forsterite. Whatever 
the reason is, it should be noted that the 40.5 /im feature is rather strong with respect to 
the 43 iim feature if it is purely due to enstatite. This implies that more materials might be 
present. 

There are some other features present around 26, 29 and 32 fim. These three posi- 
tions together with the 40.5 fim. feature could fit with the features of diopside (CaMgSi206) 
measured in the laboratory (Koike et al. 2000). The slight wavelength shift around 26 and 
40.5 yum with respect to the laboratory, and, in addition, the lack of a clear feature around 
21 /im makes the identification less conclusive. But these differences might be explained 
by temperature and grain size and shape effects as is also observed for other silicates, like 
forsterite (Bowey et al. 2002; Molster et al. 2002a,b). 

Another mineral that might explain some of the features is anorthite. Properties of 
anorthite have been measured in the laboratory (Chihara et al. 2003; Keller et al. 2003). 
The profile presented in Fig. 9 is for a sample of Angg in Chihara et al. (2003). It matches 
with the feature at 26 /xm, the long wavelength side of the 40 /xm complex and the extra 
feature just before the 33.7 //m forsterite feature. This would be the first identification of 
anorthite outside the solar system. The lack of anorthite features in the ISO-spectrum 
around 18 jim might be explained by low temperature, or interstellar absorption. 

We would require future spectra beyond 45 /im to confirm the presence of anorthite and 
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diopside, which have stronger bands in these wavelength. We should also stress here that we 
are dealing with an artificial continuum and this can have some infiuence on the strength of 
the features. 



5. Discussion 

The mid-infrared images of IRAS 16279 show an elongation in the PAH band and a 
rectangular shape in the Nil. 9 and Q-band images. At the center, the PAH-to-continuum 

ratio decreases. There is some resemblance to the Red Rectangle (HD 44179) (Waters et 
al. 1998). In the Red Rectangle, the silicates are thought to be located in a disk and the 
PAHs are in the perpendicular polar flows. The images of IRAS 16279 strengthen this link 
between morphology and mixed chemistry as proposed by Waters et al. (1996) and Molster 
et al. (1999). 

The temperature of the crystalline silicate bands in HD 44179 is about 135 K (cnstatite; 
Molster et al. (2002b)), while it is below 110 K in IRAS 16279. Nevertheless, both objects 
show a lower temperature in the silicates than in the PAH region. In fact, all crystalline 
silicates in post-AGB stars tend to show low temperatures (100-250 K) (Molster et al. 2002b). 

In the case of IRAS 16279, the model suggests that the PAH emission in the SW and 
NE regions (Fig.l) cannot be shielded by the dusty oxygen-rich regions. The oxygen-rich 
material should be configured as an torus rather than a shell, allows radiation to leak towards 
the SW and NE directions. The model explains the PAH excess at 2 arcsec: the decreasing 
temperature suppresses the continuum but not the PAH bands. 

This structure, with a low-density carbon-rich region and an obscured, dense oxygen- 
rich torus, agrees with the (Red Rectangle) oxygen-rich-disk scenario of Waters et al. and 
Molster et al. But the oxygen-rich region in IRAS 16279 is several times larger than seen 
in the Red Rectangle. A circumbinary disk is expected to be compact in order to store the 
oxygen-rich gas over a long time. It is not clear whether this has happened in IRAS 16279. 
The extended torus should disrupt faster than the disk of the Red Rectangle. Its formation 
may have happened relatively recent. 

The crystallization of amorphous silicates occurs via heating and subsequent cooling of 
the grain. This may occur slowly at low temperature in a long-term stable disk, under the 
infiuence of UV radiation (as in the Red Rectangle), or quickly in the AGB wind at very 
high mass-loss rates, through high temperature annealing (Waters et al. 1996; Sylvester et 
al. 1999). Assuming an expansion velocity of 20kms^^, the mass-loss rate of IRAS 16279 
was of order 5 x 10~'^MQyr~^, sufficient for high temperature annealing in the outfiow to 
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occur. Therefore, crystallization itself may occurred already in the AGB outflow. Part of 
the crystallized silicate is stored in the disk later, and part of it may remain in the outflow. 

It is not clear that the torus is stable long enough for long-term crystaUization. This 
may be solved by shocks addressed by Barker, Desch (2002) who found a 5 km shock is 
sufficient for anneahng comets. This velocity range might be possible in AGB or post- AGB 
wind interaction with the torus. In this case, the higher density of the torus is more likely 
to obtain a higher rate of crystallization, and this may compensate the short life of the torus 
than the disk. 

The crystalline silicate bands in IRAS 16279 are significantly weaker than in the BD 44179 

(Waters et al. 1998). The weakness of the features indicate that either they are less abun- 
dant than in HD 44179, or the temperature difference with the hotter amorphous silicates 
(responsible for the continuum) is larger in IRAS 16279. Although a difference in abundance 
can clearly not be excluded, the low temperature of the crystalline silicates seems to be 
supported by the absence of detectable features below 30 //m. 
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and F.K. is supported by NASA through the SIRTF Fellowship Program, under award Oil 
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Fig. 2. — The intensity profiles in declination, across the source center (corresponding to 
RA offset = 5.6 arcsec in Fig. 1). The peak intensity is normahzed to unity. The profiles of 
the N8.9 and Nil. 9 bands, which cover PAH features, are shown in bold: they have longer 
tails than the continuum bands (N9.8 and N10.4; thin lines). The profile of the calibration 
star (HD 169916) at N9.8 band is shown as the dashed fine. 
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Fig. 3. — The intensity ratios of the PAH bands with respect to the continuum (N9.8). The 
PAH is enhanced in the outflows. 
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Fig. 4. — The intensity ratio of the PAH bands with respect to N9.8 (pseudo-continuum) 
along declination through the source center (RA offset = 5.6 arcsec). The data is smoothed 
to an 0.4 arcsec grid. The approximate error levels are about 0.8-1.0. 
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Fig. 5. — The N-band spectra with TIMMI-2, compared to ISO/SWS spectra. The sht 
position of TIMMI-2 is RA offset ~4.0 arcsec, which crosses the south-west PAH enhanced 
region. Flux units of the TIMMI-2 spectra are in Jy per pixel. The ISO/SWS spectra are in 
Jy but are scaled by 1/3. The 11.2 /xm PAH feature is strong in the spectra located 0.9 arcsec 
South, and changes its shape from north to south. 
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Fig. 6. — The 11.2 /xm PAH feature after subtraction of the continuum (left). The continuum 
is estimated using a hnear fit between 10.5-11.0 and 11.9-12.4 //m. The PAH feature is 
strongest at the SW outflow (0.9" south). (Right) To examine the variation in the shape, 
the individual spectra are normalized to the peak intensity. We plot the spectra at 1.8" 
South in gray lines as references. Near the center, the feature is broader at the red wing 
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Fig. 7. — ISO/SWS spectra. The numbers show the identifications of crystalhne sihcate 
bands:'!' forsterite + plateau, '2' enstatite, '3' diopside, and '4' anorthite. Emission lines at 
6.63 and 34.8 //m could be due to CI. The dashed line is the two-component model described 
in the text. The dotted line represents the continuum component of the model spectra due 
to the carbon grains only. The spectra of PAH features are enlarged in the bottom. 
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Fig. 8. — The comparison between the PAH bands in the ISO/SWS spectra of HD 44179, 
IRAS 16279 and the reflection nebulosity near CD-42 11721 (top-to-bottom). The dashed 
hnes indicate where the IRAS 16279 spectrum differs most from the ISM PAH spectrum 
(as represented by CD-42 11721) and HD 44179. The spectra are binned to a resolution of 
A/AA = 200 
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Fig. 9. — The continuum-subtracted ISO/SWS spectra (solid line). The dashed line shows 
the fitting of the crystaUine bands (Molster et al. 2002b). The dust properties are measured 
by Koike et al. (1999, 2000); Chihara et al. (2003). 



